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We study the interplay between ordered and chaotic dynamics at the critical 
point of a generic first-order quantum phase transition in the interacting boson 
^ ' model of nuclei. Classical and quantum analyses reveal a distinct behaviour of 

^ ' the coexisting phases. While the dynamics in the deformed phase is robustly 

f**^ ^ regular, the spherical phase shows strongly chaotic behavior in the same energy 

intervals. The effect of collective rotations on the dynamics is investigated. 



Quantum phase transitions (QPTs) are structural changes occurring at 
' zero temperature, resulting from a variation of parameters in the quantum 

CNI ■ Hamiltonian. They have become a topic of great interest in diverse many- 

body systems, e.g., atoms, molecules and nuclei. The abrupt changes in the 
the system's ground state affect the nature of the underlying dynamics and 
can lead to the emergence of quantum chaos. In the present contribution, 
we examine this effect in generic (high-barrier) first-order QPTs in nuclei, 
focusing on the role of the barrier separating the two coexisting phases.^ 

We employ the interacting boson model (IBM),^ widely used in the 
description of quadrupole collective states in nuclei, in terms of a system of 
N monopole (s) and quadrupole (d) bosons, representing valence nucleon 
I pairs. A geometric visualization of the model is obtained by an energy 

surface, V{l3,'y), which serves as a Landau potential. The latter is defined 
by the expectation value of the Hamiltonian in the intrinsic condensate 
state \(3,r,N) = (7V!)-V2[rJ(/3, 7)]^|0), where rt(/3,7) = -i=[/3cos74 + 
sin 7-^(^2 + 0^12) + \/2 — The quadrupole shape parameters in the 

intrinsic state characterize the associated equilibrium shape. QPTs have 
been extensively studied in the IBM."^'^ To construct a critical Hamiltonian 
it is convenient to resolve it into intrinsic and collective parts, ^ 

iJcii = -ffint + Hcol ■ (1) 

The intrinsic part (Hint) determines the potential V{l3,j), while the col- 
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lective part {Hco\) contains kinetic rotational terms which do not affect 
its shape. For a first-order critical point, the two parts of the full critical 
Hamiltonian (iJcri) can be transcribed in the form® 

Hint = hPUM ■ P2{M , (2a) 



Hcol = C3 



Co (3) - 6nd 



C5 



Co(5) - 4nd 



C6 



(2b) 



Here Ps^/^o) = /9os^4 + - W is the rf-boson (total) 

number operator and Cq denotes the quadratic Casimir of the group G 
as defined in Ref. 5. Barred parameters imply scaling by N{N — 1). For 
(3o > 0, Hcii annihilates both the spherical s-condensate, \(3 = 0, 7; N), and 
the deformed condensate, |/3 = /3e = a/2/3o(1 + /3o)"^^^7 = 0;-^) which 
correspond to the two coexisting shape phases of the nucleus. 

The classical limit of the IBM Hamiltonian is obtained through the use 
of Glauber coherent states and taking A'' ooJ The resulting classical 
image; of Hcri contains complicated expressions, including square roots of 
polynomials in the coordinates /3,7 and their conjugate momenta Pi3,Pj. 
Setting all momenta to zero, leads to the potential 



V„i/h2 = + i(l - /3o')/3" - i/3o ^2^/33 cos 37 . (3) 

The potential has a spherical minimum at /3 = degenerate with a prolate- 
deformed minimum at (/3 = /3o, 7 = 0) both at energy t4iin = 0. The barrier, 
separating the two minima, is located at /3 = [1 — (1 + /3o)"^^^]^^^ ^^^d has 
a height Vb = /i2[l - (1 + PlY^'^Y 1'^- The limiting value at the domain 
boundary is K;ri(/S = V^, 7) = ^2- The potential can also be expressed in 
Cartesian coordinates x = 13 cos j and y = /Ssiwy. 

Chaotic properties of the IBM have been extensively studied,^ albeit, 
with a simplified Hamiltonian, giving rise to an extremely low, hence non- 
generic, barrier. To study the effect of the barrier, we consider the classical 
dynamics associated with ifint (2a) constrained to L = 0. The Poincare 
sections for /3o = 1-0, 1.3, 1.5, are displayed in Fig. 1. The three cases corre- 
spond to "low" , "medium" and "high" potential barriers Vb/h2 = 0.04, 0.10, 
0.16, (compared to Vb/h2 = 0.0009 in previous works'"). The sections are 
plotted at energies Ei = Vb/i, E2 = Vh and Es = 41^. In all cases, the mo- 
tion is predominantly regular at low energies and gradually turning chaotic 
as the energy increases. However, the classical dynamics evolves differently 
in the vicinity of the two wells. The family of regular trajectories near the 
deformed minimum forms a simple pattern of concentric loops around a 
single stable (elliptic) fixed point. The trajectories remain regular even at 
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Fig. 1. Poincare sections at energies (panels a^-bj^-cj^, k = 1-3) and Peres lattices of 
(A'^ = 80, L = 0) eigenstates of f^int (2a) overlayed on the classical potential Vcvi{x, y = 0) 
for /3o = 1.0, 1.3, 1.5 and /i2 = 1. Notice the well-separated regular and chaotic dynamics 
associated with the deformed {x > 0) and spherical (x = 0) minima, respectively. 

energies E ^ Vb (panels 03-63-03). In contrast, near the spherical mini- 
mum, a more complex regular dynamics at low energies (panels ai-61-Ci), 
turns chaotic at a much lower energy (E « Vh/3) and complete chaoticity 
is reached near the barrier top. The clear separation between regular and 
chaotic dynamics, associated with the two minima, persists all the way to 
E = Vb, (panels 02-62-02). In general, the regularity is more pronounced for 
higher barriers (larger Pq). These attributes are present also in the quantum 
analysis in terms of Peres lattices,® formed by the set of points {xi,Ei}, 
with Xi — y^2{i\nd\i) /N, and Ei the energy of the eigenstate The Peres 
lattices for L — eigenstates of Hint (2a) are shown in the bottom row in 
panels (a-b-c), nested within the Landau potential V{x,y — 0). For each 
/3o, one can clearly identify several regular sequences of states localized in 
and above the respective deformed wells. A close inspection reveals that 
their x^-values lie in the intervals of a;-values occupied by the regular tori 
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Fig. 2. Peres lattices for eigenstates with A'^ = 50, /3o = 1-35, h2 = 1 and L = 0,2,3,4 
of Hint (2a) (panel a) and additional collective terms (2b) involving 0(3), 0(5) and 0(6) 
rotations (panels b-c-d), respectively. Notice the well-defined rotational bands {K = 0, 
L = 0, 2, 4) and {K = 2, L = 2, 3, 4) formed by the regular states in the deformed phase 
shown in panels (a, b, c), which are distorted in panel (d). 

in the Poincare sections. In contrast, the remaining states, including those 
residing in the spherical minimum, do not show any obvious patterns and 
lead to disordered (chaotic) meshes of points at high energy E > Vh- 

Fig. 2 displays combined Peres lattices for eigenstates with N — 50, 
/3o = 1-35 and L = 0, 2, 3, 4 of iJint (2a) (panel a) supplemented with three 
different additional collective terms (2b) {ca, C5, cg} = {1, 0, 0}, {0, 1, 0} and 
{0, 0, 1} in panels (6,c,d), respectively. The regular sequences of L = eigen- 
states, mentioned previously in relation to Fig. 1, are seen to be bandhead 
states of (ivT = 0, L = 0, 2, 4) rotational bands of states with nearly equal 
values of {fid)- Similarly, sequences of L = 2, 3, 4 states form K = 2 bands. 
Such ordered band structures are not present in the chaotic parts of the 
lattice. The ordered band structure is well preserved by the 0(3) and 0(5) 
terms but is suppressed by the 0(6) term. 
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